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Abstract 

The BESS-Polar collaboration has recently performed a precise measurement 
of the local antiproton flux which is consistent with a pure secondary production 
of antiprotons. We constrain a possible primary component originating from dark 
matter pair-annihilations. We derive limits on the annihilation cross section which 
are stronger than or comparable to those from the PAMELA satellite experiment 
for dark matter masses up to 200 GeV. Especially, we exclude thermal WIMPs with 
masses in the range 3-20 GeV if they annihilate dominantly into quark pairs unless 
their cross section is velocity suppressed. 
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1 Introduction 



Antimatter searches provide a powerful tool to identify exotic sources of particles in 
our galaxy. In combination with other direct or indirect detection techniques they may 
eventually allow us to reveal the nature of dark matter. Recently, the BESS-Polar II 
experiment has measured the local antiproton flux in the range 0.2 — 3.5 GeV |T]. This 
energy regime is of utmost importance as it could contain imprints from the evapora- 
tion of primordial black holes [2] or the pair-annihilation of weakly interacting massive 
particles (WIMPs) especially with masses of (9(10 GeV). Such light WIMPs have been 
proposed as a resolution for anomalies observed by the direct detection experiments 
DAMA [2111], CoGeNT [5] and CRESST |6]. They have also been considered as the 
source of a possible gamma ray excess [7H9] and radio filament structures [10] near the 
galactic center and, recently, also in the context of the isotropic radio emission [TT]. None 
of these hints is unambiguous: the direct detection data - if interpreted in terms of dark 
matter - are in tension with the experiments CDMS [12] and XENON |13[[Ti] while 
the indirect observations may find an astrophysical explanation. Nevertheless, the light 
WIMP hypothesis requires further investigation. In the same spirit, collider data [T514T9] . 
the cosmic microwave background |20[I2T]. gamma rays and neutrino observations 

at Super-Kamiokande |27] were used to constrain the properties of light WIMPs. 

Earlier antiproton searches by BESS [281129]. AMS |30] and PAMELA |3l] have been 
analyzed in the context of light supersymmetric dark matter [32HSZ] and for the specific 
example of a 10 GeV WIMP annihilating into bottom pairs |38) . A model-independent 
study, however, seems now in order, especially as the quality of data has improved sub- 
stantially with the BESS-Polar II antarctic flight. But the high precision of the measure- 
ment also requires a profound understanding of the secondary antiproton background 
which arises from the spallation of cosmic rays on interstellar matter. 

Therefore we will recalculate the secondary antiproton flux taking into account im- 
provements in the determination of the antiproton source term and the cosmic ray propa- 
gation parameters. We will then evaluate the primary contribution from dark matter an- 
nihilation considering all relevant hadronic final states. This includes a careful treatment 
of processes like reacceleration and energy losses which substantially affect antiproton 
fluxes at low energies but have been neglected in many previous works. Finally, based 
on the BESS-Polar II antiproton search, we will present limits on the annihilation cross 
section of dark matter. 

2 Secondary Antiproton Background 

While the possibility of a primary antiproton source in our galaxy is highly speculative, 
the presence of so-called secondary antiprotons is well established. These originate from 
spallations of cosmic rays on the interstellar gas. Within error bars all antiproton searches 
performed so far are consistent with a purely secondary origin. Owing to the high pre- 
cision of the BESS-Polar II data we decided to reevaluate the secondary antiproton flux 
using the 2- zone diffusion model of Donato et al. jSH] • Compared to the original work we 
improve the calculation by including a new determination of propagation parameters |40] 
and by using the updated antiproton source term given in [41] . 



1 



2.1 Source Term 



Galactic cosmic rays which mainly consist of protons (H) and helium (He) may create 
antiprotons by inelastic interactions with the interstellar gas in the galactic disc. The 
dominant reactions are H,He + HismjHcism p + X where the index ISM stands for in- 
terstellar matter. The secondary source term which describes the differential p production 
rate per volume, time and energy takes the form 

oo 

q-{T) = 2 J2 fdT' (^) nA,3M <^'a(T') . (1) 

A=H,He ^ ^ /A+AisM^P+X 

Here $a denotes the cosmic ray fiuxefl of protons and helium which can be extracted 
from |12], riAisM their interstellar densities (tt-Hism — 0.9 cm~'^, riHejsM — ^^ism/^^ [39]). 
The differential cross section for the reaction A+Aism — ^ p + X is expressed in terms 
of the kinetic energies T' and T of the incoming nucleus and the outgoing antiproton 
respectively. The threshold energy is i?th = 6 rUp and the factor of two on the right-hand 
side accounts for p production by anti-neutron decay. 

The calculation of the source term requires a reliable parameterization of the an- 
tiproton production cross sections. The p-p interaction is well measured and an analytic 
expression for the cross section cth+h^p+x was provided in |13|H^. The interactions 
involving helium were first treated in form of a simple energy-independent enhance- 
ment factor (see e.g. |45 y 46 ] ). later a model-based evaluation using the DTUNUC Monte 
Carlo program was performed [3S1I1Z]. In this work we prefer the semi-analytic approach 
of jmHH] where cross sections are obtained by fitting the parameterization of Kali- 
novski et al. [49] to the available experimental data. We consider this method more 
reliable towards lower p energies where the DTUNUC model reaches the edge of its va- 
lidity [41| . The source term can be extracted from figure 9 in [JT], a good fit in the range 
T = 0.1 - 100 GeV is given by 

with ci = 0.98, C2 = -0.72, cg = -0.021, C4 = 0.023 and C5 = -0.0021. 



(2) 



2.2 Propagation 

The antiproton flux induced by the source term qp is determined through the diffusion 
equation. Including all processes relevant for cosmic ray propagation and assuming steady 
state the latter can be written as [50] 

V{-K VNp + V, Np) + drihot Np - Kee drNp) + r,„„ Np = qp , (3) 

where Np denotes the antiproton space-energy density. The first term on the left-hand 
side accounts for antiproton diffusion by interactions with inhomogeneities in the galac- 
tic magnetic field. Convection is included through the galactic wind velocity Vc. The 

^In principle the proton and helium fluxes change with the distance from the galactic center. We 
find, however, that the inclusion of the radial dependence only marginally affects the local antiproton 
fiux. Therefore we take $h,Ho to be homogeneous in the galactic disc. 
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last three terms on the left-hand side treat energy losses, diffusive reacceleration and 
annihilations. 

To simplify the diffusion equation we follow the approach of Donato et al. which 
we will briefly review in the foUowingj^ The halo where diffusion and convection occurs 
is approximated by a cylinder of half height L and radius R equal to that of the galactic 
disc, i.e. i? ~ 20 kpc. As a boundary condition it is imposed that the antiproton density 
vanishes at the edge of the halo. The diffusion parameter K is assumed to be homoge- 
neous over the halo, using magnetohydrodynamics considerations it can be written in 
the form [S3] 

where Kq is a normalization constant, 5 the power law index, /3 and p the antiproton ve- 
locity and momentum respectively. The galactic wind velocity Vc is taken to be constant 
and pointing away from the galactic disc. 

Energy losses, reacceleration and annihilations as well as the source term are confined 
to the galactic disc. Neglecting the thickness of the disc one has to multiply them by 
2 h 5{z) in order to keep a proper normalization. Here z parameterizes the distance from 
and h ~ 0.1 kpc the half-height of the galactic disc. The term fetot includes ionization. 
Coulomb and adiabatic energy losses as well as reacceleration. It can be taken from j5] 
The energy diffusion coefficient is also related to reacceleration, one finds 



where Va is the Alfven speed of magnetic shock waves in the galactic disc. One further has 
to consider the disappearance of antiprotons through annihilations with the interstellar 
hydrogen or helium. We use Tann = ('^Hkm + 4^^^ '^HejsM) ^ann P with 



_ J 661 mb X (l + 0.0115 (^)"° ''' - 0.948 ((^)°-°''') T < 14.6 GeV , 
[36mbx(^)"-' T> 14.6 GeV. 

Finally antiprotons may lose energy through inelastic (non-annihilating) scattering with 
the interstellar gas which leads to a redistribution of their energies. This effect can be 
included by introducing a tertiary source term in the galactic disc 

(oo 

'(7) 

where primed (unprimed) quantities refer to the antiprotons before (after) scattering 
while r denotes the radial distance from the galactic center. The non-annihilating cross 
section can be extracted from |31]. For the solution to the diffusion equation with the dis- 
cussed approximations we refer to appendix lA.li The interstellar antiproton flux follows 
from 

^f = ^PN,. (8) 



^For a full numerical solution to the diffusion equation see [?ni52| . 
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The solution Np depends on the five free parameters Kq, 6, L, Vc and Va which can partly 
be fixed by observing the nuclear composition of cosmic rays. In | ^ 1 ^ 137] the configu- 
rations which correctly reproduce the boron to carbon (B/C) ratio were determined. 

Unfortunately, the B/C ratio does not considerably constrain the size of the diffusion 
halo L. While this uncertainty only mildly affects the secondary antiproton flux it is 
important for a possible primary component which we will discuss later. The situation 
has slightly improved in the last years as there now exist several independent hints for 
L > 4 kpc arising e.g. from the observation of radioactive isotopes in cosmic rays (see 
discussion in |38| and references therein). Further, a very recent study of radio data 
shows a clear preference for L ~ 4 kpc while disfavoring smaller values of L j^S]. We 
therefore adopt L = 4 kpc in the following and take the remaining parameters from a 
new B/C analysis |10] (NORM configuration in tabled]). To be conservative we will later 
also consider a smaller L = 3 kpc where we adjust the other parameters according to 
figure 5 in |4U) (SMALL configuration in table [T]). 

model 5 A'o (kpc^ ■ Myr^^) L (kpc) K (km ■ s-^) K (km ■ s'^) 
NORM OM 0.0042 4 18/7 35^5 

SMALL 0.86 0.0031 3 18.6 30.5 

Table 1: Propagation parameters consistent with the B/C ratio. 



2.3 Solar modulation 

The diffusion equation determines the interstellar antiproton flux while experiments 
measure the antiproton flux at the top of the earth atmosphere ^J*-*^. The latter is 
affected by solar modulation. To obtain one in principle has to solve a new transport 

equation which turns out to be quite delicate as transport parameters change with time 
correlated to the solar activity. Further, the magnetic fleld of the heliosphere has a 
complicated structure and is only partly accessible to experiments. We use the force- 
fleld approximation to calculate the TOA antiproton flux [59|[60] 

(T) ~ 2mgT + IS (r + 0) . (9) 

P ^ ^ 2mp(T + 0) + (T + 0)2 ^"^^ ^ ^ 

In this simple approach particles and antiparticles are modulated in the same way and 
the value of the force fleld (p can be determined through observation of the proton flux. 
By taking the interstellar proton flux from |l2] and extracting the TOA proton flux 
from [BT], we flnd = 0.5 GV. The BESS-Polar II data where taken in a period of 
lowest solar activity where the predictions of the force-fleld approximation and more 
sophisticated models of solar modulation converge. While this gives confidence that our 
treatment is sufficient, some remarks seem in order. 

In general also drift effects play an essential role for the transport of charged par- 
ticles in the heliosphere [62 | [63 ] . The latter are charge-dependent and therefore distin- 
guish between protons and antiprotons. If one compares the antiproton measurements 
by BESS-Polar II and PAMELA [H] one finds that the fiux determined by PAMELA is 
slightly higher [T]. This difference - if it has its origin in solar modulation - can hardly 
be explained within the force field approximation as also the p/p ratio of PAMELA is 
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Figure 1: TOA antiproton flux measured by BESS-Polar II together with the predicted secondary 
flux from our calculation and from Donato et al.. Solar modulation is included through the force 
field approximation with cj) = 0.5 GV. 

larger |61) . While BESS-Polar II was operating in a very short time period around the 
solar minimum (December 2007 and January 2008) PAMELA took its data from 2006 
to 2008, i.e. partly at intermediate solar activity. Away from the solar minimum drift 
models like predict an increase of the p/p ratio which could resolve the small 

discrepancy between BESS-Polar II and PAMELA. A similar effect is also visible in ex- 
perimental data on the e/p ratio measured by the Ulysses spacecraft (see e.g. [M] and 
the discussion in [65]). 

2.4 Comparison with Experimental Results 

Having discussed the antiproton production and propagation we can now turn to a com- 
parison with experimental data and earlier work. As can be seen in figure [1] the secondary 
antiproton flux from our calculatio deviates considerably from that of Donato et al. |39| 
although our calculation is based on the same propagation model. The reason is in part 
that we took the propagation parameters from the new B/C analysis [IQ] which espe- 
cially suggests a higher galactic wind velocity Vc- This leads to a slight decrease of the 
local flux as more antiprotons can escape via convection. Further, we used a different 
source term compared to Donato et al. which is smaller especially towards low energies 
(see discussion in section l?!T]) . As shown by the BESS-Polar II collaboration (figure 3 
in p[]) the flux of Donato et al. is to shallow to well reproduce their data at low energies 
even if one allows for an arbitrary normalization factor. 

If we compare our flux with the BESS-Polar II measurement we find good consistency 
regarding the spectral shape. On average, our flux is slightly below the measured flux. 
A x^-analysis gives x^/d.o.f. =2.1 and x^/d.o.f. = 0.86 if we would normalize our flux 
by a factor of 1.1. A small underestimation of the flux can easily arise by a variety 
of systematic effects. We have neglected antiproton production on heavier nuclei which 

•^The secondary antiproton fluxes for the NORM and SMALL configurations are virtually indistin- 
guishable. Therefore, they are not discussed separately. 
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may increase the antiproton flux by ~ 5% [H]. Further, the interstellar densities of 
H, He, the propagation parameters as well as our treatment of solar modulation are 
subject to uncertainties. It is thus evident that secondary antiprotons alone are capable 
to explain the experimental data. This was also found in jl]. Nevertheless, we will take the 
conservative view that some fraction of the measured flux could still be of primary origin. 
In the following we will constrain the contribution to $p from dark matter annihilation. 



3 Antiprotons from Dark Matter Annihilation 

Dark matter annihilation can be an efficient source of antiprotons in our galaxy. We 
will now derive the primary antiproton flux trying to make only few assumptions on 
the nature of dark matter. As we restrict our analysis to the BESS-Polar II experiment 
we only consider the mass range < 200 GeV. For heavier dark matter particles the 
PAMELA satellite experiment is more sensitive and we refer the reader to |66) . 



3.1 Model-independent Approach 

Antiprotons can be generated if the annihilation products of dark matter particles x 
involve quarks or gauge bosons. We will consider 

XX dd, ss, cc, bb, WW, ZZ (10) 

and assume 100% annihilation into one channel. One can easily rescale our results with 
the corresponding branching fraction. Leptonic final states are omitted as they do not 
give rise to an appreciable antiproton production]^ 

The primary antiproton source term reads 

Here dNj^ / dT denotes the antiproton energy spectrum per annihilation for the channel 
under consideration (/ = u, d, s, c, 6, W, Z). It is determined with the PYTHIA Monte 
Carlo (version 8.1) |71] . For the dark matter density we use a Navarro- Frenk- White 
profile 

with Tq = 8.5 kpc and = 24.4 kpc. We checked that the primary antiproton fiux is 
virtually insensitive to the choice of the profile. For the local dark matter density we 
take Po = 0.39 GeVcm"'^ from |72) . a similar value was also found in [73]. Note, however, 
that this quantity is subject especially to systematic uncertainties |74] . 

The dark matter annihilation cross section (cTannf )o averaged over the current velocity 
distribution is a free parameter. In case dark matter is produced thermally as suggested 



*Note, however, that leptophihc annihilation may give rise to antiprotons through electroweak 
bremsstrahlung |67H70| . 
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Figure 2: Propagation function for different choices of the propagation parameters. The NORM 
and SMALL configurations are used in this work. They were derived from a more recent B/C 
analysis |40) compared to MIN, MED and MAX. The shaded band refers to the energy range 
of BESS-Polar IL 

by the WIMP hypothesis it is related to the dark matter density. One has to bear in 
mind that the velocity distribution of WIMPs at freeze-out is different from the current 
one. However, in a wide class of models a-^^-^v is velocity-independent and then (crann'y)o 
is fixed as [7^ 

{a^nnv), ^ 10-2^ cm^ s"i X J_ ^ , (13) 

where Tp is the freeze-out temperature. The effective number of degrees of freedom at 
freeze-out g*{Tp) can be taken from |76] . 

3.2 Primary Antiproton Flux 

The primary antiproton flux follows from the diffusion equation ([3]) with the source term 
given by (fTTj) . The main difference compared to secondary antiprotons arises from the 
fact that sources are distributed over the whole diffusion halo and not just located in 
the galactic disc. If we neglect for the moment low energy effects on the antiprotons - 
namely energy losses, reacceleration and tertiaries - the primary antiproton space energy 
density can be written as (see appendix IA.2I) 

iV^"^C(r0,^)7^o, (14) 

where TZq denotes the propagation function. Note that this formula only provides an 
estimate, in our analysis we will always use the full solution to the diffusion equation 
which we discuss in appendix IA.2I Nevertheless, the function TZq nicely illustrates the 
dependence of A'^"™ on the propagation parameters. It is shown for the NORM and the 
SMALL configurations in figure [2l To allow for comparison with earlier work we also 
depict TZq for the commonly used MIN, MED and MAX configurations [32] which were 
derived from an earlier B/C analysis |57]1^ 

^The propagation function for the MIN, MED and MAX configurations was extracted from |77| . 
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Figure 3: Comparison of secondary (sec) and primary (prim) antiproton fluxes for the case 
of a dark matter particle annihilating into bottom quarks with mass = 8 GeV (left) and 
= 30 GeV (right). A cross section of ((Tann^^)o = 3 x 10~^^ cm^ s~-^ is assumed. Also shown 
are the primary fluxes without taking into account low energy effects (prim approx). 

At high energies cosmic ray propagation is dominated by diffusion and the differences 
between the five configurations mainly stem from the choice of the halo size L. The two 
configurations used in this paper (NORM and SMALL) approach the MED configuration 
as they have the same or a similar value of L. However, NORM and SMALL decrease 
more rapidly towards low energies. In this regime the process of convection rapidly gains 
importance, and - as suggested by the new B/C analysis |40| - NORM and SMALL 
assume a larger convective wind than MIN, MED and MAX. In the energy range of 
BESS-Polar II we will therefore obtain primary antiproton fluxes in our analysis which 
are similar as in the MIN configuration. 

In figure in] we depict the primary antiproton fluxes (NORM configuration) expected 
from a thermal WIMP annihilating into bottom quarks for two different masses. To illus- 
trate the importance of energy losses, reacceleration and tertiaries we show the primary 
flux with and without taking into account these low energy effects. For a first impression 
of the sensitivity of BESS-Polar II, we have also included the secondary background and 
the data points. It can be seen that the 8 GeV WIMP is clearly inconsistent with the 
measurement, while the 30 GeV WIMP is still viable. 

4 Limits from BESS-Polar II 
4.1 Our Analysis 

As discussed in section 12.41 the secondary antiproton flux alone is capable to explain 
the BESS-Polar II measurement, there is no need to include a primary component. We 
therefore do not try to improve the fit to the data by invoking dark matter annihilations. 
We rather want to determine the maximal rate of annihilations still consistent with 
the experiment. For this we will use a modified x^-metric which only punishes a given 
configuration if the total (primary + secondary) antiproton flux overshoots the data. We 
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Figure 4: Limits on the dark matter annihilation cross section for the NORM and SMALL 
configurations. The shaded band refers to the annihilation cross section of a thermal WIMP in 
the case of s-wave annihilation (see text). 



define 



Xmod 



El ^ <J>TOA ^ ^data 

) ^ '^P,^ > ^P,^ ' (15) 

I " ^TOA ^ ^data 

P,i P,i ■ 



Here denotes the predicted flux averaged over the bin i, ^p'f'^ the measured flux 

and (Tj the statistical error in the bin. Note that Xmod does not follow an ordinary 
probability distribution. For 29 d.o.f. the 95% upper limit on the primary flux corresponds 
to xLd = 25.5. 



4.2 Results 

In flgurelHwe provide 95% upper limits on the dark matter annihilation cross section for 
all hadronic channels. We show the results separately for the two propagation parameter 
sets of table [U The NORM conflguration assumes a diffusion halo size L = 4kpc compared 
to L = 3 kpc in the SMALL conflguration, the latter resulting in more conservative 
constraints. We also depict the annihilation cross section of a thermal WIMP for the 
case of s-wave annihilation. The latter is obtained from (1131) where we take m^/Tp to be 
in the range 20 — 25 and include an additional 20% uncertainty on the cross section if 
the freeze-out happens around the QCD phase transition. 

It can be seen that a thermal WIMP with s-wave annihilation into quarks is excluded 
in the mass range 3 — 20 GeV even for the SMALL conflguration. At smaller masses the 
antiproton production threshold is approached and predictions become uncertain. The 
limits for the bottom and charm channels are weaker than those for light quarks especially 
towards low WIMP masses. The reason is that the heavy quarks induce fewer antiprotons 
and the spectrum is more compressed towards lower energies, i.e. more antiprotons reside 
below the energy threshold. 
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Constraints also arise for heavier dark matter particles. The primary antiproton flux 
from a 100 GeV WIMP annihilating into W pairs is e.g. still to ~ 30% in the energy range 
of BESS-Polar II (0.2 — 3.5 GeV). Nevertheless, the limits become weaker as the number 
density of dark matter scales inversely with its mass. Therefore, above 30 GeV only 
scenarios with a non-thermal dark matter annihilation cross section can be constrained. 

5 Conclusion 

We have studied the production and propagation of antiprotons in the light of new exper- 
imental data by BESS-Polar II. After having recalculated the cosmic ray induced flux we 
conflrmed that the measurement is compatible with a purely secondary origin of antipro- 
tons. We determined the primary flux which would arise from dark matter annihilations 
considering all relevant annihilation channels and taking into account updates on the 
propagation parameters. The latter especially show a trend towards higher convective 
winds which considerably reduces the primary antiproton flux. Nevertheless, because of 
the high precision of the BESS-Polar II data, we were able to provide strong constraints 
on dark matter annihilations. We found that a thermal WIMP with = 3 — 20 GeV 
is excluded if it dominantly annihilates into quarks unless the cross section is velocity 
suppressed. While this statement relies on rather conservative astro- and nuclear physics 
assumptions it is impossible to exclude all sources of uncertainty. In this light it is im- 
portant to note that our result is conflrmed independently by observations of gamma 
rays from dwarf galaxies Pl] and of the cosmic microwave background (SDIEI]- As the 
uncertainties affecting the three analyses are widely uncorrelated the exclusion of light 
WIMPs with hadronic annihilations seems relatively robust. 

This implies that the signals of the direct detection experiments DAMA, CoGeNT 
and CRESST - if interpreted in terms of dark matter - require non-hadronic channels, 
a non-thermal cross section or p-wave annihilation. The flrst of the named options is 
again subject to strong constraints especially from the Super-Kamiokande neutrino tele- 
scope |27] . 

Although the BESS-Polar II experiment took its data at low energies it can still 
provide important limits on dark matter particles with masses up to (9(100 GeV) as 
antiprotons typically only carry a small fraction of the energy of the mother particle. 
Thermal WIMPs with > 30 GeV are still viable, but certain models which involve 
boost factors for dark matter annihilation or non-thermal cross sections can be con- 
strained. 

A Solving the Diffusion Equation 
A.l Secondaries 

To get rid of the radial part of the diffusion equation one can Bessel expand the space- 
energy density and the source terms 

oo oo 

Np{r,T) = Y,Np,iT)Jo {(-) , g,-(r, T) = Y.qpAT)Jo (O.^) • (16) 

1=1 i=l 
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Here i? ~ 20 kpc denotes the radius of the galactic disc and Q the i-th zero of the Bessel 
function Jq. The Bessel coefficients for the secondary source term are determined by 

R 

r \ 

(17) 







R 



As discussed in section [2m we assume that g|*''^(r, T) = g^'^'^(T) 6{R — r) where 6 denotes 
the Heaviside function and g|'^'^(T) is given by ([T]). 

Plugging ( !T6|) into (|3]) and performing the approximations described in section [2l2] the 
axial part of the diffusion equation can be solved analytically. One arrives at a differential 
equation in energy which reads (at z = 0) 



Ap,{Np, - iV°J + 2hdT{btot Np, - Kee drNp,) = 2hq^ 



with 



and 



Ap, = 2/ir,„„ + K + KS, coth (^^^ , Si=^^ + A^. (19) 
= ■ (20) 

Pi* 

The equation for the antiproton space energy density has to be solved numerically. In 
the high energy regime (T ^ 10 GeV), energy losses, reacceleration as well as tertiaries 
can be neglected and Np^i approaches Np^^. 



A. 2 Primaries 



For the case of primary antiprotons, the diffusion equation can be solved in the same 
way as for secondaries, the only difference being the distribution of sources. It turns out 
that the differential equation 0181) is still valid, however, iV? • has to be defined differently. 
One finds (see e.g. [32]) 

^^°. = C"^(r0,T)7^, (21) 

with 




Here we have Bessel expanded the dark matter density distribution p^ir, z) and denoted 
the Bessel coefficients by (p^)j. The latter are determined analogous to the q^pf in (fT7|l . 

Again the equation for Np^i has to be solved numerically. In the high energy limit 
Np^i approaches Np^ and the local antiproton space-energy density reads 

oo 

Np ~ g?™(r0, T) 7^0 , 7^0 = 7^, Jo (o^) • (23) 

i=l 

Note, however, that for antiproton energies T < 10 GeV this equation should only be 
used as a first estimate. In our analysis we will consider the full numerical solution to 
the diffusion equation and only use (123|) for illustrative purposes. 
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